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A search for a diffuse flux of astrophysical muon neutrinos, using data collected by the ANTARES neutrino 
telescope from December 2007 to December 2009 is presented. A (0.83 x 2ii) sr sky was monitored for a total of 
334 days of equivalent live time. The searched signal corresponds to an excess of events, produced by astrophysical 
sources, over the expected atmospheric neutrino background without any particular assumption on the source 
direction. Since the number of detected events is compatible with the number of expected background events, a 
90% c.l. upper limit on the diffuse i/^ flux with a spectrum is set at E'^^<jQ% = 5.3 x 10~* GeV cm"'^ s~^ sr~^ 
in the energy range 20 TeV - 2.5 PeV. Other signal models with different energy shape were also tested and some 
rejected. 



1. Introduction 

The ANTARES high-energy neutrino telescope 
is a three-dimensional array of photomultiplier 
tubes (PMT) distributed over 12 lines installed 
deep in the Mediterranean Sea, each line includ- 
ing 75 PMTs [T]. A neutrino telescope in the 
Northern hemisphere includes the Galactic Cen- 
tre in its field of view and is complementary to 
the IceCube Antarctic telescope [2]. 

The main goal of the experiment is the 
search for high-energy neutrinos from astrophys- 
ical sources. If the sensitivity of point source 
search techniques is too small to detect neutrino 
fluxes from individual sources, it is possible that 
many sources could produce an excess of events 
over the expected atmospheric neutrino back- 
ground. In this proceeding the search for very- 
high energy extraterrestrial muon neutrinos from 
unresolved sources is presented using data col- 
lected by the ANTARES telescope from Decem- 
ber 2007 to December 2009. 

Atmospheric muons and neutrinos are the main 
sources of background in a neutrino telescope. 
The former can be suppressed by applying re- 
quirements on the direction of the events, the 
latter is an irreducible background. As the spec- 
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trum of cosmic neutrinos is expected to be harder 
than that of atmospheric neutrinos, the signal we 
are looking for corresponds to an excess of high 
energy events in the measured energy spectrum 
without any particular assumption on the source 
direction. 

Electrons (in the so-called "leptonic models") 
or protons and nuclei ( "hadronic models" ) can be 
accelerated in astrophysical processes. Hadronic 
models [3. predict that the energy produced in the 
sources is carried away by cosmic rays, 7-rays and 
neutrinos. A benchmark flux for the measure- 
ment of diffuse neutrinos is the Waxman-Bahcall 
(W&B) upper bound [1]. Using the CR observa- 
tions at EcB^ ~ 1019 eV {EIj^^cr. ~ 10"* GeV 
cm~^s~isr~i) the diffuse flux of muon neutrinos 
is constrained at the value: 

El^^ < 4.5/2 X 10"* GeV cm"^ sr^^ s"^ (1) 

(the factor 1/2 is added to take into account neu- 
trino oscillations). 

2. Neutrino tracking and energy recon- 
struction 

Muon neutrinos are detected via charged cur- 
rent interactions: Vfj^ + N ^ fj, + X . The arrival 
times and the amplitudes of the Cherenkov light 
signals detected by the PMTs [S] are used to re- 
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construct the trajectory of muon neutrinos and 
to estimate their energy. 

The track reconstruction algorithm defined in 
[B] is based on a hkehhood fit that uses a detailed 
parametrization of the probability density func- 
tion for the photon arrival times. The track posi- 
tion and direction, the information on the number 
of hits {Nhit) used for the reconstruction and a 
quality parameter A are the main outputs of the 
algorithm. A is determined from the likelihood 
and the number of compatible solutions found by 
the algorithm itself. A can be used to reject badly 
reconstructed events. For > 10 TeV, an an- 
gular resolution for muon neutrinos better than 
0.3° is accomplished by the ANTARES detector. 

2.1. Monte Carlo simulations 

The Monte Carlo (MC) simulation tools [715] 
include the production of Cherenkov light, the 
generation of the optical background caused by 
radioactive isotopes and bioluminescence present 
in sea water, and the digitization of the PMT sig- 
nals. In particular, the PMT simulation also in- 
cludes the probability of a detected hit giving rise 
to an afterpulse. The simulation of afterpulses is 
critical when the energy estimator defined in §2.21 
is applied to MC events. The afterpulse probabil- 
ity was measured in laboratory using ANTARES 
[H] and NEMO [10, PMTs and it was confirmed 
with deep-sea data. Upgoing muon neutrinos and 
downgoing atmospheric muons have been simu- 
lated and stored in the same format used for the 
data. 

Signal and atmospheric neutrinos. MC 

muon neutrino events have been generated in the 
energy range 10 < E'^ < 10^ GeV and zenith 
angle between 0° < 9 < 90° (upgoing events). 
The same MC sample can be differently weighted 
to reproduce the "conventional" atmospheric neu- 
trino spectrum (Bartol), oc E^^ '^ at high energies 
[TT] . and the expected astrophysical signal spec- 
trum, oc E~^. The normalization of the signal 
flux is irrelevant when defining cuts, optimizing 
procedures and calculating the sensitivity. Here 
a diffuse flux test signal is deflned equal to: 

El $^ = 1.0 X lO"'^ GeV cm-^ s"^ sr"^ (2) 

Above 10 TeV, the semi-leptonic decay of 
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Figure 1. Definition of the variable Ri on the i-ih 
PMT. In this example, both ARSO and ARSl are 
fired; after the integration dead-time, both chips col- 
lected light again. In this example, Ri = 4. 

short-lived charmed particles D — > K + ^ + 
becomes a significant source of atmospheric 
"prompt leptons" . The Recombination Quark 
Parton Model (RQPM) is used in this simula- 
tion, since it gives the largest prompt contribu- 
tion among the models considered in [12] . 

Atmospheric muons. The ANTARES trig- 
ger rate pLS^ is dominated by atmospheric muons 
that represent the main background for a neu- 
trino telescope. A small fraction (approxi- 
mately 5%) of triggered downgoing muons is mis- 
reconstructed as upgoing; their rejection is a cru- 
cial point in this analysis. 

The MUPAGE package [M] was used to sim- 
ulate atmospheric muon samples. One year of 
equivalent live time with a total energy Ex ^ 
1 TeV and bundle multiplicity m = 1 1000 
was generated. The total energy Et is the sum 
of the energy of the individual muons in an at- 
mospheric muon bundle. Triggered ANTARES 
events mainly consist of multiple muons originat- 
ing in the same primary CR interaction [15] . 

2.2. Energy dependent variable 

The only way to separate atmospheric and as- 
trophysical neutrinos is through a discrimination 
based on the energy. An original energy estimator 
is defined, which is based on hit repetitions in the 
PMTs due to the different arrival time of direct 
and delayed photons (Fig. [1]). Direct photons 
are emitted at the Cherenkov angle and arrive 
at the PMTs without being scattered. Radiative 
processes contribute to energy losses linearly with 
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Figure 2. Reconstructed muon energy using tlie R 
energy estimator as a function of the true muon en- 
ergy. The displayed events are selected after the qual- 
ity cuts defined in tj3.1l Only neutrino- induced muons 
are selected. 

the muon energy for E'^ > 1 TeV. The resulting 
electromagnetic showers produce additional light. 
Photons originating from secondary electromag- 
netic showers or scattered Cherenkov radiation 
arrive on the PMTs delayed with respect to the 
direct photons, with arrival time differences up to 
hundreds of ns [3]. The fraction of delayed pho- 
tons increases with the muon energy. 

The signal produced by the PMTs is processed 
by two Analogue Ring Sampler (ARS) [TBJ which 
digitize the time and the amplitude of the sig- 
nal (the hit). They are operated in a token ring 
scheme. If the signal crosses a preset threshold, 
typically 0.3 photo-electrons, the first ARS inte- 
grates the pulse within a window of 40 ns. If 
triggered, the second chip provides a second hit 
with a further integration window of 40 ns. After 
digitization, each chip has a dead time of typi- 
cally 250 ns. After this dead time, a third and 
fourth hit can also be present. 

The number of repetitions Ri for the i-th PMT 
is defined as the number of hits in the same PMT 
within 500 ns from the earliest hit selected by the 
reconstruction algorithm (Fig. [T]). In most cases, 
Ri —1 or 2, but it could be also 3 or 4. The mean 
number of repetitions in the event is defined as 

R — ^ ^' , where Npmt is the number of PMTs 

NPMT ' 

in which hits selected by the tracking algorithm 
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Number of expected events for MC and data. 
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Data 
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80 
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24 


1.0 ■ 10"* 
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116 


20 





Expected events in 334 days of equivalent live time 
for the three MC samples: atmospheric muons, 
atmospheric neutrinos (Bartol-|-RQPM) , astrophys- 
ical signal from eq. [2j and data. Reco: at the 
reconstruction level; Upgoing: reconstructed as 
upgoing; f^'-level: after the first-level cuts; 2"'*-level: 
after the second-level cut. 

are present. For a muon neutrino sample, R is 
linearly correlated with the log of the true muon 
energy Etme in the range from 10 TeV (i? ~ 1.26) 
to 1 PeV (i? ~ 1.73). R can be used to estimate 
the muon energy Ereco, see Fig. [51 The distribu- 
tion of \og{ErecolEtrue) has & FWHM=0.8. The 
resolution is comparable or better with respect to 
other energy reconstruction algorithm |17) . 

This energy estimator is robust because it does 
not depend on the number of active PMTs and 
on non-linear effects on charge integration. 

3. Cosmic neutrino signal selection 

The data collected from December 2007 to De- 
cember 2009 are analyzed. In this period, the de- 
tector configuration changed several times with 9, 
10 and 12 active lines. For this reason, three dif- 
ferent detector configurations, based on the num- 
ber of active lines and optical modules, were re- 
produced in MC simulations. Data runs are se- 
lected according to the data-quality requirements 
explained in [S]. The total live time is 334 days: 
70 days with 12 lines, 128 days with 10 lines and 
136 days with 9 lines. 

3.1. Rejection of atmospheric muons 

As described in ^2.1[ the rejection of atmo- 
spheric muons is a crucial point in the search for 
a cosmic neutrino signal. This contamination can 
be strongly suppressed by applying requirements 
on the geometry of the events and on the track re- 
construction quality parameter A. Two different 
levels of cuts are defined in order to remove the 
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Figure 3. Scatter plots of the reconstruction quality parameter A versus the number of hits Nhit for atmospheric 
muon (left) and signal neutrino (right) MC simulations after the first-level cuts. The pink line represents the 
second-level cut described by eq. [S] The color code is in units year~^. 



contamination of mis-reconstructed atmospheric 
muons from the final sample. 
First-level cuts. Events are selected according 
to these criteria: (i) upgoing particles with recon- 
structed zenith angle Orec < 80° (corresponding 
to 0.83 x27r sr); (m) A > -6; (Hi) Nhit > 60; (iv) 
reconstruction with at least two lines. The first- 
level cuts reduce the rate of mis-reconstructed 
events by almost 3 orders of magnitude, as in- 
dicated in Table [TJ 

Second- level cut. The remaining atmospheric 
muons have the quality parameter A which on av- 
erage decreases with increasing Nhu- Fig. |3](left) 
shows the correlation between A and Nhu for at- 
mospheric muons. In order to completely remove 
the expected rate of mis-reconstructed events in 
the MC sample, a cut value A* is defined as a 
function of Nhit- 

r -4.59 - 5.88 • lO-^'Nh.t N^t < 172 
1 -5.60 NMt > 172 ^ 

Removing all events with A < A*, the atmo- 
spheric muons are completely suppressed. Inde- 
pendent MC atmospheric muon simulations us- 
ing CORSIKA (see details in [S]) confirm that 
the maximum contamination in the final sam- 
ple is less than 1 event /year. As can be seen 
in Fig. |3] (right), the signal is highly preserved 
from the second-level cut. The effects of the first- 
and second-level cuts on signal and atmospheric 
neutrinos are also given in Table [TJ 



3.2. Discrimination from atmospheric neu- 
trinos 

A cut on the energy dependent variable R, de- 
fined in §2.21 is used to separate the diffuse fiux 
signal from the atmospheric background. The 
optimal cut value is obtained through a blinding 
procedure on MC events, without using informa- 
tions from the data. The numbers of expected 
signal (us) and background {nt) events are com- 
puted as a function of R. Then, calculating the 
so-called Model Rejection Factor (MRF) defined 
in [TB], the best cut is obtained and used as the 
discriminator between low energy events, domi- 
nated by the atmospheric neutrinos, and high en- 
ergy events, where the signal could exceed the 
background. After the optimization of all the 
parameters, the observed data events {uobs) are 
revealed (un-blinding procedure) and compared 
with the expected background for the selected re- 
gion of R. If data are compatible with the back- 
ground, the upper limit for the signal flux is cal- 
culated using the Feldman- Cousins method [19] 
at a 90% confidence level (c.L). 

The cumulative distributions of the R vari- 
able for atmospheric neutrino background (Bar- 
tol-l-RQPM) and diffuse flux signal (eq. [2]) are 
computed for the three conflgurations of the 
ANTARES detector and the corresponding live 
times. The MRF is calculated as a function of R 
using these cumulative distributions. The mini- 
mum found for R — 1.31 determines the cut value 
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Figure 4. Distribution of the R parameter for the 
134 neutrino candidates in the 334 days of equivalent 
live time. Points represent data, the filled histogram 
is the atmospheric neutrino MC (Bartol model only). 
The dashed line represents the maximum contribu- 
tion (RQPM) of "prompt" neutrinos. The MC pre- 
dictions are not normalized to the data. The signal 
normalized at the upper limit (eq. |4} is shown as a 
full line. The cut at R = 1.31 is indicated as a vertical 
line. 

for the energy dependent variable ^2U^. Assum- 
ing the Bartol (Bartol+RQPM) atmospheric 
flux, 8.7 (10.7) background events and 10.8 sig- 
nal events are expected for R > 1.31. The central 
90% of the signal is found in the neutrino energy 
range 20 TeV < < 2.5 PeV. 

4. Upgoing neutrino candidates 

4.1. Low energy events R < 1.31 

Events surviving the second-level cut are up- 
going neutrino candidates. The first step of the 
un-blinding is to reveal the events with R < 1.31. 
In this region, 125 events are found. A compar- 
ison with MC predictions is shown in Fig. |4] as 
a function of R. The events with R > 1.31 in 
Fig. [4] are uncovered after the final un-blinding 
of the data sample. MC predictions are lower 
by ~ 20% with respect to the detected events. 
Bartol atmospheric neutrino MC predicts 104.0 
events with i? < 1.31, and Bartol + RQPM pre- 
dicts 105.2 events. The discrepancies between 
predicted and measured events are well within the 
systematic uncertainties of the absolute neutrino 



flux at these energies (25-30%) [H]. 

The number of expected background events 
with R > 1.31 is 8.7 for Bartol MC only Most 
prompt models described in jl2J give negligi- 
ble contributions; the RQPM model predicts the 
largest contribution of 2.0 additional events with 
respect to the conventional Bartol flux. An aver- 
age over all the considered models gives a contri- 
bution of 0.3 events. A combined model of Bartol 
flux plus the average contribution from prompt 
models is adjusted with the data/MC normal- 
ization factor obtained in the R < 1.31 region. 
Hence the number of expected background events 
for R > 1.31 is 10.7. 

4.2. High energy events and upper limit 

The MC simulations have been tested and com- 
pared with data. In particular, the R distribu- 
tions show a reasonable agreement both for atmo- 
spheric muons 20^ and for atmospheric neutrinos 
in the low energy region R <1.31 (c.f. Fig. 2]). 
As a consequence, the signal region with R > 1.31 
was un-blinded and 9 high-energy neutrino can- 
didates are found. 

Systematic uncertainties on the expected num- 
ber of background events in the R > 1.31 region 
are evaluated considering: (i) the contribution 
of prompt neutrinos, estimated as events; 
(m) the uncertainties from the conventional neu- 
trino flux, that depend mainly on the uncertainty 
on the absolute flux as a function of the energy 
and on the spectral index, evaluated to be ±1.1 
events. The uncertainties on the detector effi- 
ciency (angular acceptance of the optical module, 
water absorption, scattering length, trigger sim- 
ulation and the effect of afterpulses) amount to 
5%: they affect the detection both of signal and 
background neutrinos in the high energy region. 

The number of observed events is compatible 
with the number of expected background events. 
The 90% c.l. upper limit on the number of sig- 
nal events /igo%(no&s, «b) for Uobs = 9 observed 
events and nt — 10.7 ± 2 background events in- 
cluding the systematic uncertainties is computed 
with the method of [21]: the value fJ-go%{nb) = 5.7 
is obtained. The upper limit on the diffuse flux 
is given by $90% = <I>^ • ^90%/ns: 

E^%o% = 5.3 X 10"^ GeV cm^^ sr'K (4) 
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Table 2 



Tested flux models. 


Model 


R* 


Nmod 


^-^90% 


tJ-90% / ^mod 








(PeV) 




MPR 


1.43 


3.0 


0.1^10 


0.4 


P96p7 


1.43 


6.0 


0.2-MO 


0.2 


S05 


1.45 


1.3 


0.3-^ 5 


1.2 


SeSi 


1.48 


2.7 


0.3^20 


0.6 


Mpp + p7 


1.48 


0.24 


0.8^50 


6.8 



Astrophysical flux models, the value of the R* which 
minimizes the MRF, the expected number of events 
Nmod, the energy range Ai?9o% in which the 90% of 
events are expected, and the ratio fJ.9o%/Nmod- See 
[20| and references therein. 



Figure 5. The ANTARES 90% c.l. upper limit for a 
diffuse z/^ + flux obtained with this analysis, 
compared with the results obtained by other experi- 
ments and theoretical predictions. See [3] and refer- 
ences therein. 

This limit holds for the energy range between 20 
TeV to 2.5 PeV. The result is compared with other 
measured flux upper limits and theoretical predic- 
tions in Fig. [5] 

Some theoretical predictions of cosmic neutrino 
fluxes with a spectral shape different from 
are also tested. For each model a cut value R* is 
optimized following the procedure in ^'d.2\ Table 
[2] gives the results for the models tested: the val- 
ues of R* , the numbers Nmod of signal events 
for R > R*, the energy intervals where 90% of 
the signal is expected, the ratios between Hgo% 
(computed according to [19]) and Nmod- A value 
of f^9o%/ Nmod < 1 indicates that the theoretical 
model is inconsistent with the experimental result 
at the 90% c.l. [20]. 

5. Conclusions 

Using data from 334 days of equivalent live 
time collected with the ANTARES telescope, a 
search for a diffuse flux of high energy cosmic 
muon neutrinos was made. The 90% c.l. upper 
limit on the diffuse i/^ flux with a E~'^ spectrum is 
set at E'^^go% = 5.3 x IQ-^ GeV cm'^ g-i sr^^ 
in the energy range 20 TeV - 2.5 PeV. Other sig- 
nal models with different energy shape are also 
tested and some of them excluded at the 90% c.l.. 
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